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Abstract. Basing on the analogy between the coherent states of light and 
separable states of N bosons, we demonstrate that the violation Cauchy-Schwarz 
inequality for any-order correlation function signals the entanglement among the 
constituent particles. Rather than restricting to the correlations between the 
positions of particles, we consider the broadest set of measurements allowed by 
quantum mechanics. Our result is general - it applies to any quantum system 
of bosons, even when the number of particles is not fixed, provided that there 
is no coherence between different number states. We also demonstrate that the 
compact expression for the separable state of bosons can be used to relate some 
known metrological quantities to the particle entanglement in a very simple way. 
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Although the foundations of quantum and classical physics are much different, it is 
often difficult to construct a simple criterion of non-classicality of a particular system. 
A prominent example of a quantum phenomenon is the ability of particles to exist 
in superpositions of states. The most well known manifestation of such superposition 
is the Young double-slit experiment for massive particles, which confirms their wave 
character and the ability to coherently interfere. A more subtle example of genuinely 
quantum phenomenon is the Einstein-Podolsky-Rosen (EPR) paradox [T], which was 
introduced using a pair of entangled particles, separated far apart. A consequence of 
quantum mechanics is that a measurement of some physical quantity of one of the 
particles affects the other at an instant. This has no classical counterpart, as shown 
by John Bell [21 El 11]. Nowadays, the Bell inequalities, tested in many experiments 
laiellTlIHlIllllollIIllIllllsllIlllIS], are analyzed in the context of the entanglement 

[13- 

Systems of entangled particles [m nH] are potentially vialable to quantum 
information [19], teleportation [201 El] or ultra-precise metrology [22l [23]. The 
definition of an entangled state is discriminative. Namely, it is based on the definition 
of the separable (non-entangled) states as being a mixture of product states of N 
particles [211 EH 1^; be., 

Qsep = ® ® ■ ( 1 ) 


Here pfs are non-negative weights that add up to unity. Entangled are all those 
states which cannot be represented in the form O. A consequence of this indirect 
definition is that it is usually not straightforward to construct a criterion or a measure 
of particle entanglement. Usually, one must find some physically measurable quantity 
A., which for separable states is bounded (from above or below) by a critical value 
^ 0 - If some state gives the value of which breaks this bound, then ^ is a criterion 
for particle entanglement. Trivially, such criterion detects entanglement of only those 
states, which break the bound, therefore most ,/l’s are not universal. 

Quantum interferometry provides a good example of such criterion. Namely, 
consider a two-mode state g oi N particles passing through an interferometer. During 
the propagation, a phase 6 is imprinted between the modes. If the phase is estimated 
in a series of m ^ 1 experiments, the precision of the estimation is limited by the 
Cramer-Rao lower bound [23 EH]: 


A6» ^ 


1 1 


( 2 ) 


Here, Fq is the quantum Fisher information, which depends on g and the type of 
applied interferometer [23- For all separable two-mode states Fq ^ N [23: therefore 
A0 is limited by the shot-noise. Nevertheless, there exist some “usefully” entangled 
states, for which Fq exceeds this critical value. Thus Fq = N plays the role of 
mentioned above. 

For atomic interferometers, the useful particle entanglement has been achieved 
by means of two-body interactions present in ultra-cold systems. Usually, such 
correlations are associated with the spin-squeezing of a two-mode sample [301 E3 
[32l 123 EH E3 E3 E3 Ell E3 • Alternatively, in a process which resembles the down- 
conversion, the interactions drive the scattering of pairs of entangled atoms from a 
coherent source, such as a Bose-Einstein condensate (BEG) [40l |4T1 l42l 1^ |44] . 
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In this work, we first show in Section[2]that the separable states of bosons formally 
resemble the classical states of electromagnetic field from Eq. ©• In Section [3] we 
derive the general expression for the second order correlation function of separable 
states. Rather then focusing on correlations between the position or momentum 
measurements, we consider the most general measurement operators allowed by 
quantum mechanics. In Section |4] we argue that the violation of the Cauchy-Schwarz 
inequality (CSI) for the second order correlation function can be considered as a 
simple and useful criterion for entanglement of bosons. We derive the CSI for 
generalized measurements and then focus on the correlations between the position 
measurements (Section 14.11) . We show, how this criterion can be extended to higher 
order correlations ('Section l4.2l) and systems, where the number of particles is not fixed 
f Section |4.31) . These Sections are an extension of our previous work [4^. In Section [S] 
we further exploit the compact expression for the separable bosonic states to simplify 
the derivations of some well known interferometric criteria for particle entanglement. 

2. Separable state of N bosons 

In this Section, we relate the coherent states of light to separable states of N 
bosons. Generally speaking, the non-classical electromagnetic field is that consisting 
of individual photons. This statement can be quantified if the density matrix of the 
electromagnetic field is expressed using the P-representation [liiiT] 

g = yp$|$)(4>|P(4>). (3) 

Here, |<i)) is a (possibly multi-mode) coherent state of light, defined by the relation 

f(+)(r)|$) = $(r)|$). (4) 

The operator f(+)(r) implies the positive-frequency part of the electromagnetic field 
£ (r), while the symbol 2?$ in Eq. ([3]) is the integration measure over the set of complex 
fields 4>. In general, the state of light is called classical if the measurement outcomes 
can be explained in terms of classical electromagnetic waves. It happens when the 
P-representation can be interpreted as a probability distribution, which occurs when 
V is normalized and 

yT>$P($)F($) > 0 (5) 

V 

holds for any volume V and all non-negative functions E(<I>). In the opposite case, the 
field is not a mixture of coherent states and we call it quantum. 

For a system of N particles (distinguishable or not), a particle-separable state is 
defined by Eq. ([IJ. Any hermitian operator can be diagonalized, hence the A^-body 
density matrix reads 

Q = ^'PkHk){'>fk\- ( 6 ) 

k 

Here, fc is a composite index constructed from the index i from Eq. o and the indeces 
labeling the eigenstates of all N one-body matrices. Similarly, the probabilities Vk 
are products of pfs and the eigen-values of . Each \'il)k) is an Wbody pure and 
separable state, which reads 

Hk) = \(l>k^)®---®\4>^k'^)- 


(7) 
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The upper index labels the particles and implies that the single-body eigen-vectors 
\4>l^'^) and |(/)^" might differ for nf^n'. 

We now restrict these general considerations to bosons. In such case the states 
I'i/'fc) must be symmetrized with respect to the particle interchange. However, after 
the symmetrization of |'0fc), the particles become entangled and consequently g is not 
separable anymore. Clearly, the only pure, symmetric and separable state must be 
a product of N identical one-body states. Therefore, the upper index n is no longer 
necessary and the state from Eq. © transforms into 

\i^u) = \<t>kr^ = \f^cN). ( 8 ) 

Therefore, the set of density matrices of separable bosons is formed by mixtures of 
states (15|) and reads 

g = Y.Vk\cfk\N){cl>k-,N\ (9) 

k 

and when the labeling is done by a continuous parameter (f 

g = (10) 

Here, symbol Vcj) denotes the measure of integration over the set of states 
The probability distribution V{4>) is normalized and for all positive-defined functions 
^ 0 a relation 

^ 0 ( 11 ) 

V 

holds for any volume V [J. The Equations and (HH closely resemble Equations 
m and ([n|) for electromagnetic fields, but while the former pair is constructed with 
coherent states of light, which posses coherences between different particle-number 
states, the latter is valid for fixed N. If the density matrix cannot be written in the 
form cni) with P{(j)) satisfying Eq. (fTTll . the state is not separable so it is particle- 
entangled. In the following Sections we exploit this knowledge to derive some criteria 
for the separability of bosons. 

3. Second order correlation function 

At this point, we formulate a criterion for the separability of a system of N bosons, 
based on their second-order correlation function. We define this function using the 
positive-operator valued measures (POVMs) denoted by These objects form the 
broadest set of measurement operators allowed by quantum mechanics [BUT]. They 
are non-negative and form a complete basis, i.e. 

= i ( 12 ) 

where f labels the detection outcomes. Most commonly, the measurements are 
modeled by the projection operators |0(^|, which form a subset of POVMs, but for the 
moment we keep the discussion fully general. We also assume that the measurements 
do not correlate particles, which means that the POVMs are single-body operators. 

I General form of separable states of N bosons has been already analyzed in different context, see 
for instance Ha HU. 
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The second order correlation function is the rate of the coincidence of the results 
f and defined as 




N 




Y, Tr gMOEjia 


(13) 


The labels i and j indicate that the operators Ei{f,) and Ejif) act on different 
particles, while the summation ensures that all coincidences contribute to the 
correlation function. For indistinguishable particles these operators do not depend 
on the labels of the particles, therefore the indeces i and j can be dropped. The sum 
gives a combinatory prefactor a 2 = N{N — 1) and upon using Eq. (|T(I1) we obtain 

= a2 (w) 

where {E{£)) = {(j)\ E{ff)\(j)). This is a general expression for the two-body correlation 
function for separable states of bosons. Any criterion for particle entanglement, which 
is based on the has to check if the second order correlation function can be written 
in the form of Eq. m- 


4. Cauchy-Schwarz inequality as an entanglement criterion 


First, we derive a simple criterion for the particle entanglement using the CSI for 
the [Sni EU EH EH EH EH Ezj- Consider two subsets Xa and of the 
measurement outcomes and the second order correlation functions integrated over a 
pair of such regions 

= JdfJdf a, (15a) 

Gif = jdfJd^’G^^H^,a, ( 156 ) 

Xb Xb 

Qff = jdfjdf'G^^\f,0- (15c) 

Xb 

The integrals (IlSal) and (11561) describe the effective local correlations within the sets 
Xa and Xb respectively. The last line (fT^ describes the cross-correlations between 
the two sets. Note that for separable states, for which the G*-^^ is given by Eq. da, 
we obtain 

g\f = Jv^M^)fM^ ( 16 ) 

where i and j are either a or b and 

/,/,(</.) = Jdfm))- ( 17 ) 

Xi/, 


Since the functions fa/bi4>) are non-negative, the CSI for integrals reads 


E><f> fb{(l>) < 


77</./2(</.) 




(18) 
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Therefore, we arrive at an inequality 


C 2 


Gas 


^ 1 , 


^( 2 )^( 2 ) 

yao 


(19) 


which is satisfied for all separable states of N boson^. Any state breaking this bound 
must be particle-entangled. 


4-1. Special case - position measurements 


We focus on an experimentally relevant case, when the POVMs represent the 
measurement of position. For the sake of clarity, we label the measurement outcomes 
with r and r' rather then f and f. The POVM representing the idealized position 
measurement (with perfect accuracy and efficiency of the detector) is a projection onto 
the position state, i.e. E{r) = |r)(r|, which clearly satisfies conditions (1121) . The second 
order correlation function G'^^^(r, r') measures the coincidence rate for detecting one 
particle at position r and the other at r'. The correlation functions in Eqs (I15al) - (ll5c|) 
are integrated over a and 6 , which in this case denote two regions of space. 

According to the CSI (Eq. (fT^ l. for separable states of N bosons, the cross¬ 
correlation between these regions cannot exceed the geometric average of the local 
correlations. On the other hand, the enhanced cross-correlation might be signaled 
by the reduced fluctuations of the difference of the particle count between the two 
regions - the so-called number squeezing. To link the CSI and the number squeezing 
we express the second order correlation function using the bosonic field operator ^(r), 

G(2)(r,r') = {^4'+(r)4'+(r')^(r')^'(r)J> . (20) 

The atom number operator in each region is defined as 

na/b= ydr^'f(r)4'(r). (21) 

Xa/h 


The fluctuations of the atom number difference operator n = fia — fit, normalized to 
the shot-noise level ntot = (na) + {ns) give the number-squeezing parameter 
((An) 2 ) 




( 22 ) 


^tot 

Using Equations (12(1 . (1^ and the integrated correlation functions (Eqs (I15all - (ll5cp l 
we obtain 


rj^ = 1 + 


^( 2 ) 

zfaa 


q(2) 

^bb 


-2Gif-(n}^ 


^tot 


(23) 


The system is number squeezed if 77 ^ < 1, which sets a constraint for the second part 
of the above equation, i.e. 

^"<1 ^ Gil^+Gl,^,^-2gif-(n}^<0. (24) 

Although this condition resembles the CSI (IT9l) . the presence of the term (n)^ turns 
out to be crucial. To picture this, consider a separable pure state of N particles 


e = \4'o;N){(j)o;N\ 


(25) 


§ The volumes of the two regions in Eqs. (I15 qII - II15cII can be chosen arbitrarily small, and in fact the 
CSI also applies to correlations between two points. However, from the experimental point of view 
it is usually advantageous to increase the signal by accumulating the data from substantial volumes. 
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divided into two unequal parts a and b. For this coherent state the integrated second 
order correlation functions are related to the atom counts in the two regions, i.e. 


Sia = (na)^(l - jf) 


C(2) 

^bb 


= 1 - 


= - jf) and = {na){nb){l - 

{{fia) - {nb)f 


Therefore 


Ar2 


< 1 


(26) 


although the state is separable and does not violate the CSI. This example shows 
that the number-squeezing can be present in systems of bosons even in the absence of 
particle entanglement. 

However, when there is symmetry between regions a and b such that (h) = 0 and 
G^aa = G^bb obtain that 

= 1 - 2 ~ . (27) 

^tot 

In this case 


rf <l ^ C 2 > 1, 


(28) 


so in such symmetric cases the number squeezing and CSI are strictly related. 
The violation of the CSI implies the number squeezing which signifies the particle 
entanglement. 


^.2. Extension to higher order correlations 

The CSI discussed so far is based on the second order correlation function. Here we 
demonstrate that the arguments leading to Eq. m can be extended to higher order 
correlations. Consider an 2TO-th order correlation function between the measurement 
outcomes .. .^ 2 m- In analogy to expression for the from Eq. (fldll we obtain 
that for separable states 

.. .,6m) = «2m jvefVicf) (E(6)) . . . (E(6m)), (29) 

where a 2 m = ™ second order correlation function, we 

introduce two subsets of the outcomes 6 integrate m variables of Eq. (HH) over Xa 
and other m over Xb obtaining 

(30) 

Similarly, GimT^^ and gH"^^ are the 2m-th order correlation functions integrated 2m- 
times over a or b, respectively. The CSI applied to Eq. o gives 

Q{2m) 

C2m = , ^ 1 

/^{2m) ^(2m) 

yGka 'Gib 

for any separable state of bosons. 


( 31 ) 







Criteria for particle entanglement in many-body systems of bosons 


4-2.1. Application to twin-Fock state We argue that the strength of the CSI criterion 
increases when m grows. To this end, consider a twin-Fock state of N particles in two 
modes a and b 


1 


m 




(32) 


This state is particle entangled, due to the indistinguishability of the constituent 
bosons. By identifying the modes with the two regions Xa and Xf,, we analyze the CSI 
from Eq. m as a function of m for the twin-Fock state. The local correlation 

function is 




(f)! 




(33) 


and symmetrically for • On the other hand, the cross correlation function is 


-f{2m) _ 


= ((at6t)™(aS)™) = 


(f)! 




We substitute these results into Eq. (EB and obtain 

(f-2m)! 


— 


(34) 


(35) 


For given N, this is a growing function of the correlation order 2m. In Fig. [I] we 
plot C 2 m from Eq. (I55|) . N = 100, 250, 500 and 1000 and the correlation order 
2m = 2,4,6,8. For small values of the parameter e = 2m/N, Eq. (1351) can be 
approximated by 

C2m - exp(e2Ar/2), (36) 


which is in excellent agreement with the results shown in Fig.[I] Clearly, the deviation 
from the classical limit C 2 m = 1 is higher when the ratio of 2m to N increases, i.e. when 
the amount of information about the system extracted from the high-order correlation 
function is large. This result might be of direct experimental relevance, since already 
a 6-th order correlation function among thermal atoms was measured [58] . 


4-3. Extension to systems with fluctuating N 

We will now demonstrate that the Cauchy-Schwarz criterion for particle entanglement 
can be successfully applied to systems where the number of particles differs from 
shot to shot. To account for these fluctuations we impose the super-selection rule, 
which excludes coherences between states with different numbers of particles [59] . 
Consequently, a general separable state of bosons is a mixture 

OO 

where each qn is given by Eq. OT, while the statistical weights pn add up to unity. 
The second order correlation function calculated with the density matrix ((571) is a 
simple generalization of Eq. (HI and reads 

OO 

,c') = X! 

N^O 



{N)jvcfVNmm){HO) (38) 
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2 4 6 8 

Correlation order 2 m 


Figure 1. (Color online) The coefficient C 2 m ^ function of the correlation order 
for different numbers of particles in the twin-Fock state. Clearly, the deviation 
from the classical limit C 2 m = 1 (shown here by a horizontal black line) increases 
with growing m. 


In this expression, pN, a 2 {N) and Vn{4i) depends on N. On the other hand, the 
average values of the single-particle operators E calculated with the one-particle states 
\(j)) do not depend on the number of particles. Therefore, it is reasonable to introduce 

CXD 

( 39 ) 

N=0 

This is an averaged V function, weighted with the coefficients a 2 {N). One can 
introduce two regions a and b and the integrated second order correlation functions as 
in Equations (I15all - (ll5cl) . Those functions are given by Eq. (fT6|) with a2P{4>) replaced 
with P{(l)). If VNif)) for each N satisfies condition ([5]), so will 'P{(j)), therefore, the 
CSI will be satisfied also by all separable states with fluctuating number of particles 
given by Eq. dSzl). Therefore, the violation of the CSI signals particle entanglement, 
just as when N is fixed. 

However, there is one subtle difference between these two cases. Namely, Vn 
can be partially negative for some iV, which means that in this sector particles are 
entangled. Nevertheless the averaged 'P{4>) can still be positive-defined, simply because 
the separable sectors of the density matrix might overshadow the entangled ones. Thus 
the CSI for systems with fluctuating number of particles can overlook the particle 
entanglement present in some fV-sector. The above considerations apply also to higher 
order correlation functions and the related CSI formulated in Section |421 

5. Metrological entanglement criteria 

Quantum metrology provides some useful criteria for entanglement of particles in 
two-mode systems. These criteria arise from general considerations involving the 
upper bounds for the precision of the estimation of an unknown phase 6, imprinted 
on a system upon passing through a two-mode interferometer. We show that the 
general expression for the density matrix of a collection of bosons from Eq. TO vastly 
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simplifies in the two-mode case. Using this expression, we reformulate the well-known 
proofs that the quantum Fisher information and the spin-squeezing parameter are 
criteria for particle entanglement. 

5.1. Two-mode separable states of bosons 

The two-mode states of bosons are described using the second quantization 
annihilation operators a and b. We can define the two-mode equivalent of the separable 
pure states defined in Eq. ®. We put all N bosons in the same quantum state which 
in the most general form reads 



(40) 


Here z S [0,1] measures the population imbalance between the two modes, while 
(fi S [—7r,7r] is the relative phase. The states |z, are called the coherent spin 

states. The relations 


d\z, (p; N) = \/iV-v/ze®‘^|z, ip;N — 1) 


(41a) 
(41 &) 


S|z, (/?; N) = — z\z, (p; N — 1) 


show the analogy with the states |(/>; N). Furthermore, the general two-mode state of 
bosons (in analogy to Eq. (nni) is a mixture of different states (HOI) , i.e. 


1 TT 



(42) 


0 —TT 


Naturally, for g to be particle-separable, V{z,p) must posses the properties of the 
probability distribution as in Eq. ®. With Eq. (14211 at hand, we proceed to the 
metrological criteria for particle entanglement. 

5.2. Quantum Fisher information for bosons 

When a collection of particles passes through an interferometer, a relative phase 0 is 
imprinted between the two modes. The interferometer can be usually represented 
by an evolution operator U = where h is the generator of the unitary 

transformation. Subsequently, some measurement is performed at the output to 
deduce the value of the parameter. Usually, the sequence of measurements is repeated 
TO ^ 1 times, to obtain the average value of the estimated parameter and the 
associated fluctuations A0. If the estimator is unbiased, i.e. it has the desired property 
that this average tends to the true value of the parameter, then according to the 
Cramer-Rao lower bound [27l [28] the fluctuations are constrained by 


1 1 


(43) 


A6» ^ 


Fq [p] 


Here Fq is the quantum Eisher information (QFI), which depends on both the density 
matrix g and the interferometric transformation [29) . The expression for Fq is in 
general rather complicated but for pure states g it reduces to 


FQ[p]=4((Ah)"), 


(44) 
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where Ah = h — {h). The average values are calculated in the input state g. When h 
is a one-body operator, so the interferometer does not correlate the particles, then 


FQ[g]i^N 


(45) 


for all separable two-mode states of N particles. Therefore Fq > N signals the 
entanglement between the particles. The original proof of this statement dealt with 
general separable two-mode states which can be written in the form of Eq. (HD m- 
We show how this derivation simplifies for two-mode states of bosons. Note that 
the two-mode interferometric transformations are generated by the angular momentum 
operators 


ix = ^ (d'^6 -I- 

4 = l(atS-Sta) 


(46 a) 
(466) 
(46 c) 


which form the closed Lie algebra of the angular momentum. The generic 
interferometric transformation IJ = where = n • J is a scalar product 

of a unit vector and J = {Jx, Jy, JzY^ ■ Note that Jn for any h can be obtained by 
rotating the operator, namely 




^ — J 


(47) 


where va is some vector oriented in a direction depending on n. In the evaluation of 

Eq. (H^ . these two external rotations can be moved so that they transform the 

density matrix rather than the operator J^. Since these operations do not entangle 
particles, the transformed density matrix of a separable state will still have the form 
of Eq. (l42l) . This means, that it is sufficient to show that Fq ^ N for the Jz 
transformation and any separable two-mode state. According to the above argument, 
this inequality will then hold for any 

One important property of the QFI is its convexity, namely 


Fq 




< '^P^FQ[gi]. 
i 


This property, applied to states from Eq. dH gives 


1 TT 

Fq [e] ^ jdzjdip'P{z,ip)FQ [\z,(p;N) {z,(p;N\]. 

0 —TT 


(48) 


(49) 


Now the QFI within the integral is calculated for a the pure state go = 
\z,ip;N) (z,(p;N\ and thus it is simply given by four times the variance of the Jz 
operator. With help of the relations (141 al) and (14161) we arrive at 

Fq[0o] =4((Ai,)') =4fVz(l-z) ^iV (50) 

for z G [0,1]. Using the property from Eq. we get 

Fq [g] < N (51) 

for all separable states. We conclude that Fq > TV is a criterion for particle 
entanglement. 
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5.2.1. Fluctuating N The above arguments can be extended to situations when the 
number of particles fluctuates in the two-mode system m- In such case in the absence 
of coherence between states with different numbers of particles the separable two-mode 
state of bosons is given by Eq. dSZl), where pat is given by Eq. dUD. For each sector 
of fixed TV, we can repeat the derivation leading to Eq. (ED) and using the convexity 
property (|48ll . we obtain 

OO 

(52) 

N=0 

Thus Fq > (N) singals the entanglement between the particles. 


5.3. Spin-squeezing for bosons 

The QFI is a powerful criterion for particle entanglement, which is useful for 
metrological purposes. However, it has one major drawback - it is hard to determine 
in the experiment. Quantum interferometry provides much simpler entanglement 
criteria, which do not require the full knowledge of the density matrix g. Among 
these is the spin-squeezing parameter [301IM] 


es=N 


W + W' 


(53) 


This parameter is bounded for two-mode separable states by ^ 1. Therefore 
< 1 implies the particle entanglement. As in the previous Section, we can show 
that the general proof of this fact presented in [24] simplifies for separable states of 
indistinguishable bosons. 

Note that with help of relations (141 oD and (141 &I) we obtain that 


z)cosvj 

(54a) 

z) sin (p 

(54b) 


((AJ,)") = ^ + iV(iv _ 1) y dzP{z) i 


- y dzP{z) 

We now apply the Cauchy-Schwarz inequality for integrals 


1 


^ 2 


^ N^J dzP{z)z{l ~ ^) J cos^ip 

{JyY ^ dzP{z)z{l — z) y d(pP{(f) sin^t^. 

Therefore, the denominator of Eq. (1531) can be bounded as follows 
{Xf + {Jyf J dzP(z)z(l - z). 

The spin-squeezing parameter is thus not smaller then 






N f dzP(z)z(l-z)’ 


(54c) 

(55 a) 
(55b) 

(56) 


( 57 ) 
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where the nominator is given by Eq. (154 el) . Our aim is to show that 1 or 

equivalently 


((AJ^)^) > TVy dzP{z)z{l - z). 
Using (154 dl . the above inequality gives 

I dzP{z) - 0 - J dzP{z) 



(58) 


(59) 


which is always true, since it is a variance of the variable z — With this we have 
proven that ^ I for all separable two-mode states of bosons. 


5.3.1. Fluctuating N When the number of particles is not fixed, the moments of 
the angular momentum operators are calculated using the two mode state dSZl). As a 
result, the Equations (154 all - (154 c|) are transformed into 


iX) 

{Jy) 


(N) J dz j dipP{z,(p)\/z{l - z) cos ip 
(N) J dz J dpP{z, p)^/z{l - z) sinyi 


(60 a) 
(606) 


= -^ + {{N^) - {Nf) J dzPiz) (^z - 0 
-{N)^ jdzP{z)l^z-^^ . 

The Cauchy-Schwarz inequality can again be applied to Equations (160 al) and 
giving 


(60 c) 

(I60l 


{J.? + {Jy? ^ {Nf J dzP{z)z{l - z). 


(61) 


Upon performing manipulations similar to those leading to Eq. (1591) but with N 
replaced with (N) and using {N'^) > (IV)^ obtain that 


{N) 






w 


> 1 


(62) 


for all separable two-mode states of bosons with the fluctuating number of particles. 
Therefore we proved that the properly defined spin-squeezing parameter is also an 
entanglement criterion in the non fixed-A case. 


6. Conclusions 

We argued that ample experience of quantum optics, and in particular its methods 
to tackle the notion of classical, semiclassical and quantum, provide tools to describe 
quantum correlations in many body systems. We considered a general POVM - a 
measurement operator allowed by quantum mechanics - and focused on the correlation 
functions of the outcomes of such measurement. We proved that when system is in 
a separable state, the correlations between the measurement outcomes satisfy the 
Cauchy-Schwarz inequality, and hence we drew the conclusion that when it is violated 
it can be treated as a criterion for non-classicality of bosonic states. Our results are 
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general, they apply to any collection of bosons, even when their number fluctuates from 
shot to shot. The only constraint is that the super selection rule applies, forbidding 
coherences between different number states. Finally we take advantage of the simple 
form of the separable state of bosons to express known interferometric criteria for 
particle entanglement, such as the QFI or the spins-squeezing in a simple manner. 
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